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Abstract 

Three-dimensional Ising model with the plaquette-type (next-nearest-neighbor and four-spin) 
interactions is investigated numerically. This extended Ising model, the so-called gonihedric model, 
was introduced by Savvidy and Wegner as a discretized version of the interacting (closed) surfaces 
without surface tension. The gonihedric model is notorious for its slow relaxation to the thermal 
equilibrium (glassy behavior), which deteriorate the efficiency of the Monte Carlo sampling. We 
employ the transfer-matrix (TM) method, implementing Novotny's idea, which enables us to treat 
arbitrary number of spins for one TM slice even in three dimensions. This arbitrariness admits 
systematic finite-size-scaling analyses. Accepting the extended parameter space by Cirillo and co- 
worker, we analyzed the (multi) criticality of the gonihedric model for N < 13. Thereby, we found 
that, as first noted by Cirillo and co-worker analytically (cluster- variation method), the data are 
well described by the multi-critical (crossover) scaling theory. That is, the previously reported 
nonstandard criticality for the gonihedric model is reconciled with a crossover exponent and the 
ordinary three-dimensional-Ising universality class. We estimate the crossover exponent and the 
correlation-length critical exponent at the multi-critical point as = 0.6(2) and z> = 0.45(15), 
respectively. 
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I. INTRODUCTION 



Study on surfaces spans a wide variety of subjects ranging from biochemistry to high- 
energy physics , leading to a very active area of research. In particular, the problem 
of interacting surface gas Q, Q, [s , 0] is of fundamental significance. The Savvidy-Wegner 
(gonihedric) model , |lO|, describes the interacting closed surfaces without surface 

tension. The surfaces are discretized in such a way that they are embedded in the three- 
dimensional cubic lattice, and the surface faces consist of plaquettes. The gonihedric model 
was introduced as a lattice-regularized version of the string field theory |12]. However, 
recent developments dwell on the case of three dimensions, aiming a potential applicability 
to microemulsions. 

The gonihedric model admits a familiar representation in terms of the Ising-spin variables 
{Si} through the duality transformation; namely, the plaquette surfaces are regarded as the 
magnetic-domain interfaces. To be specific, the Hamiltonian is given by the following form; 

H = Ji ^ SiSj + J2 ^ SiSj + J3 ^ SiSjSkSi, (1) 

with finely tuned coupling constants, Ji = —2k, J2 = k,/2 and J3 = —(1 — /t)/2. The 
Ising spins Si = ±1 are placed at the cubic-lattice points in three dimensions, and the 
summations ^(jj), and X][jjfc«] '^^^^ possible nearest-neighbor pairs, next- 

nearest-neighbor (plaquette diagonal) spins, and round-a-plaquette spins, respectively. The 
interfacial energy E of the gonihedric model is given by the formula E = n2 + Ann^, where 
n2 is the number of links where two plaquettes meet at a right angle (folded-link length) and 
n4 is the number of links where four plaquettes meet at right angles (self-intersection-link 
length). Namely, the surfaces are subjected to a bending elasticity with a fixed strength, 
and the the self-avoidance is controlled by the parameter k. We notice that the interfacial 
energy lacks the surface-tension term. 

Because of the absence of the surface tension, thermally activated undulations should be 

Hmoted significantly. Such feature might be reflected by the phase diagram; see Fig. [T] (a) 
ll^ . We notice that a phase transition occurs at a considerably low temperature quite 
reminiscent of that of the two-dimensional Ising model. Moreover, for large n, the phase 
transition becomes a continuous one, whose criticality has been arousing much attention: 
By means of the Monte Carlo method, Johnston and Malmini obtained the critical 



exponents v = 1.2(1), 7 = 1.60(2) and P = 0.12(1) for the self-avoidance k = 1. (Here, 
we quoted one typical set of exponents among those reported in the literature by various 
means.) The authors claimed that the exponents bear remembrance to those of the two- 
dimensional Ising model, namely, z/ = 1, 7 = 7/4 and (3 = 1/8. On the other hand, 
with the Monte Carlo method, Baig and co-worker obtained u = 0.44(2) (k = 1) and 
7/// = 2.1(1) {k = 0.5, 1). By means of the low-temperature series expansion, Pietig and 
Wegner |17l obtained a = 0.62(3), (3 = 0.040(2) and 7 = 1.7(2) (k = IV With use of the 
cluster- variation method with the aid of the Pade approximation 19, 20, Cirillo and 
co-worker obtained the estimates /? = 0.062(3) and 7 = 1.41(2). 

Meanwhile, a subtlety of the Monte Carlo simulation inherent to the gonihedric model 



was noted by Hellmann and co-worker 



22| . According to them, the relaxation to the thermal 



equilibrium is extremely slow, and such slow relaxation smears out the singularity of the 
phase transition. In order to cope with such slow relaxation (long auto-correlation length), 
they employed the histogram Monte Carlo method. However, the singularity of the phase 
transition could not be resolved satisfactorily. (See also Ref. [2i| for an alternative evidence 
of strong metastabilities.) As a matter of fact, the gonihedric model at k = reduces to 
the so-called ferromagnetic p-spin model, and the model has been studied extensively as a 



possible lattice realization of super-cooled liquids and glassy behaviors 
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28l |. In this sense, an alternative simulation scheme other than the Monte Carlo method 



is desirable in order to surmount the slow-relaxation problem and determine the critical 
exponents reliably. 

In this paper, we develop a transfer-matrix formalism, implementing Novotny's idea 
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33, 



34l |. which enables us to treat arbitrary number of spins N for one transfer-matrix 



slice. This arbitrariness admits systematic finite-size-scaling analyses. (In addition to this 
advantage, the transfer-matrix calculation yields the correlation length ^ directly. Because ^ 
has a fixed scaling dimension, succeeding effort at the finite-size-scaling analyses is reduced 
to a considerable extent.) We also accept the idea of Cirillo and co-worker who extended 
the parameter space of the gonihedric model ((H) to, 

1 — K 

Jl = -1, J2 = -j, and J3 = — . (2) 

4k, 

(Note that for j = —0.25, the parameter space reduces to that of the aforementioned orig- 
inal gonihedric model.) With respect to this extended parameter space, Cirillo and co- 
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worker claimed that the above mentioned pecuhar criticahty could be identified with a mere 
end-point singularity (multi-criticality) [2^ of an ordinary critical line of the three- 
dimensional- Ising universality class; see the critical branch of the phase diagram shown in 
Fig. n(b). Thereby, they obtained the crossover critical exponent cj) = 1-1(1) by means of 
the cluster- variation method Our transfer-matrix simulation supports their idea that 
the numerical data are well described by the multi-critical (crossover) scaling theory. We 
estimate the crossover exponent and the correlation-length critical exponent as = 0.6(2) 
and z> = 0.45(15), respectively; hereafter, we place a dot over the critical indices at the 
multi-critical point (j = —0.25). 

The rest of this paper is organized as follows. In Sec. IHl we set up a transfer-matrix 
formalism for the gonihedric model based on Novotny's idea. In Sec. IIIIl we present the 
numerical results. Taking the advantage of the Novotny formalism, we carry out systematic 
finite-size-scaling analyses. In the last section, we present summary and discussions. 

II. EXTENSION OF THE NOVOTNY METHOD TO THE PLAQUETTE-TYPE 
INTERACTIONS 

In this section, we present methodological details of our numerical simulation for the 

nnn 

gonihedric model (0). We employed Novotny's improved version j31|,|32,|3^ of the transfer- 



matrix method. This technique allows us to construct the transfer matrices containing 
arbitrary number of spins in one transfer-matrix slice; note that in the conventional 
scheme, the available system sizes N are limited for high spatial dimensions d > 3 severely. 
Actually, Novotny constructed the transfer matrices of the Ising model for d < 7 fairly 
systematically j34|. Such arbitrariness of admits systematic finite-size-scaling analyses. 

In the following, we adopt Novotny's idea to study the gonihedric model ([Q). For that 
purpose, we extend his idea so as to incorporate plaquette-based interactions. We restrict 
ourselves to the case of three dimensions d = 3 relevant to our concern. (The original idea 
of Novotny is formulated systematically for general dimensions, taking the advantage that 
only the bond-based (nearest neighbor) interaction is involved.) 

We decompose the transfer matrix into the following three components. 
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where the symbol denotes the Hadamard (element by element) matrix multiplication. 
Note that the multiplication of the local Boltzmann weights should give rise to the total 
Boltzmann factor. The decomposed parts, T^^^a) ^ j^ipianar) ^ j^(rung)^ p g^^nd for 

the Boltzmann weights for intra-leg plaquettes, intra-planar plaquettes, and rung plaquettes, 
respectively; see Fig. |2] as well. 

First, let us consider the contribution of T^^^a\ The matrix elements are given by the 
formula, 

J-ij - \I'\^\JI - S{i,l)S{i,2) S{i,2)S{i,2.) ■ ■ ■ S(i,N)S(i.,l) ' V^j 

where the indices i and j specify the spin configurations for both sides of the transfer-matrix 
slice. More specifically, we consider spins for a transfer-matrix slice, and the index i 
specifies a spin configuration {S{i, 1), S{i, 2), . . . , S{i, N)} arranged along the leg; see Fig. 
121 The factor Wg'^^^ denotes the local Boltzmann weight for a plaquette with corner spins 
{Si, . . . , 5*4}. Explicitly, it is given by the following form, 

^Ifst = (-^ ( ^(^1^2 + ^2^4 + 5-4^3 + ^3^1) + y (^1^4 + ^2^3) + y^i 52^3^4) ) ■ 

(5) 

(The denominators of the coupling constants are intended to avoid double counting.) Here, 
the parameter T denotes the temperature. It is to be noted that the component T*^'*^^\ with 
the other components ignored, leads the transfer-matrix for the two-dimensional gonihedric 
model. The other components of 7'(p''*"-'*'~) and should introduce the "inter leg" 

interactions so as to raise the dimensionality to d = 3. 

Second, we consider the component for the intra-planar interaction. It is constructed by 
the following formula, 

^(planar) ^ ^^|^pvW|,^^ (g) 

where the matrix P denotes the translation operator; namely, with the operation, a spin ar- 
rangement {S{i,m)} is shifted to {S{i,m + 1)}; the periodic boundary condition is imposed. 
An explicit representation of P is given afterward. Because of the insertion of P^ , the pla- 
quette interaction A bridges the V^-th- nearest-neighbor pairs, and so, it brings about the 
desired inter-leg interactions. This is an essential idea of Novotny's work. Crucial point is 
that the operation P^ is still meaningful, even though the power is an irrational num- 



ber 
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32 
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This rather remarkable fact renders freedom that one can choose arbitrary 



number of spins. 



An explicit representation of is given as follows [31|, |32, |33[. As is well known, 
the eigenvalues {pk} of P belong to the roots of unity like exp(i0fc) with (pk = 2iTk/N 
{k = 0,1, . . . , N — 1). The complete set of the corresponding eigenvectors are constructed by 
the formula = A^^^^ '^h=iPk^'^\^) ■ Here, the set {]$)} consists of such bases independent 
with respect to the translation operations, and A^^^ is a normalization factor. Provided that 
the eigenstates are at hand, one arrives at an explicit representation of P^; 

(^I^1j) = E(^I'^^M(<^>;c|j)- (7) 

Finally, we consider the component of T^^'^'^^\ This component is also constructed sim- 
ilarly. This time, however, we need two operations of P^, because concerns both 
sectors of i and j (both sides of the transfer- matrix slice); see Fig. |21 The elements are 
given by, 

Tlr'^ = HA ® {J\)B{{P^\^)) ® (p^Ij))) , (8) 
where the operator B acts on the direct-product space; 

N 

{{^\ 01) 10) = n (9) 



m=l 



Putting the components, T^'^f), T^p^^'"''') and T^™"^'), into Eq. ©, we obtain the complete 
form of the transfer matrix. Actual numerical diagonalizations are performed in the following 
section. 



III. NUMERICAL RESULTS 

In this section, we survey the criticality of the gonihedric model (0) for the extended 
parameter space by means of the transfer-matrix method developed in the preceding 
section. In particular, we investigate the critical branch with an emphasis on the end-point 
singularity at j = —0.25. We neglect a possible deviation of the multi-critical point from 
j = —0.25 as pointed out by the cluster- variation-method study Such deviation is so 
slight that it would not affect the multi-critical analyses very seriously We treated the 
system sizes up to A^ = 13. The system sizes A^ are restricted to odd numbers, for which 



the transfer-matrix elements consist of real numbers 
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A. Survey of the critical branch with the Roomany-Wyld approximative beta 
function 



To begin with, we survey the criticahty of the second-order phase boundary in Fig. ^(b). 
For that purpose, we calculated the Roomany-Wyld approximative beta function (3^^{T). 
We stress that the availability of /5^^(T) is one of major advantages of the transfer-matrix 
method. The Roomany-Wyld beta function is given by the following formula 3^, 

I _ ln(gjv(T)/gjv.2(T)) 

/?«^(T) = .^Q. 



«iv(T)Cjv-2{T) 

Here, C,n{T) denotes the correlation length for the system size N, The correlation length is 
readily calculated by means of the transfer-matrix method. That is, using the largest and 
next-largest eigenvalues, namely, Ai and A2, of the transfer matrix, we obtain the correlation 
length ^ = l/ln(Ai/A2) immediately. 

In Fig. El we plotted the beta function /?^^(T) for various j with the fixed self- avoidance 
parameter k = 2. The zero point (fixed point) of the beta function (3^^ {T) indicates the 
location of the critical point Tc. In Inset of Fig. IHl we plotted the phase-transition point 
Tc{j). This phase boundary corresponds to the critical branch of the phase diagram shown 
in Fig. n(b); the other phase boundaries are of first order, and the determination of them 
is out of the scope of the present P§^(T) approach. 

The slope of the beta function at T = Tc yields an estimate for the inverse of the 
correlation- length critical exponent l/u. In Fig. El we also presented a slope (dashed 
line) corresponding to the three-dimensional-Ising universality class u = 0.6294 for a 
comparison. We see that the criticality is maintained to be the three-dimensional-Ising 
universality class for a wide range of j. More specifically, for j = —0.05, 0.1, 0.25, 0.4, 0.55, 
and 0.7, we obtained the correlation-length critical exponent as z/ = 0.634, 0.641, 0.643, 
0.643, 0.642, and 0.642, respectively. From this observation, we estimated the exponent along 
the critical branch as z/ = 0.638(5) fairly in good agreement with the three-dimensional-Ising 
universality class. 

It is to be noted that, as mentioned in Introduction, at j = —0.25, very peculiar critical 



exponents have been reported so far 



Id 



19l |. The above simulation result suggests 



that such peculiar criticality should be realized only at j = —0.25 (critical end-point). This 
idea was first claimed by Ref. with the cluster- variation method. In fact, on closer 



inspection, the beta function in Fig. El shows a crossover behavior such that the slope in the 
off-critical regime is enhanced; see the regime of T — Tc > 3 at j = —0.05 in particular. It 
appears that such regime of enhancement is pronounced as j —0.25. Eventually, right at 
j = —0.25, a new universality accompanying small z>(< u) may emerge. In the succeeding 
subsections, we provide further support to this issue. 

For the region in close vicinity to the critical end-point, for instance, —0.25 < j < —0.2, 
we found that the beta function acquires unsystematic finite-size corrections; even the zero 
point of (T) disappears. In this sense, we suspect that a direct simulation at j = —0.25 
would not be very efficient. Rather, performing simulations for a wide range of j, we are 
able to extract informations concerning the end-point singularity fairly reliably. 



B. End-point singularity of the critical amplitude of ^ 

In the above, we found that the universality class of the critical branch is maintained to 
be that of the three-dimensional Ising model. A notable feature is that a crossover to a new 
universality class emerges as j — > —0.25 . In this subsection, we study this multi-criticality 
in terms of the theory of the crossover critical phenomenon. We read off the crossover 
exponent (p from the end-point singularity of the amplitude ^| of the correlation length. 
Namely, the correlation length should diverge in the form, 

^ « iV^|T-Tc|"^ (11) 

with the amplitude, 

oc A^-^+'^y'l'. (12) 

Here, the variable A denotes the distance from the multi-critical point A = j + 0.25. (It 
is to be noted that the critical point Tc depends on A as demonstrated in Inset of Fig. 01) 
The above formula is a straightforward consequence of the multi-critical (crossover) scaling 

nn 

hypothesis j29j, |30|| ; 

^ ^ \T -T^\-^X{A/\T -T^l"^). (13) 

As noted in the previous subsection, the dotted critical index stands for that right at the 
multi-critical point. 

To begin with, we determine the critical amplitude A^^. In Fig. |3J we plotted the scaled 
correlation length (T — T^jL^^^-^lT — T^]'^ for k = 2 and j = 0.3. The symbols, +, x, *, 
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□ , and ■, denote the system sizes of = 5, 7, 9, 11, and 13, respectively. The hnear 
dimension of the system L is given by L = \fN . In the plot, we postulated the three- 
dimensional- Ising universality class v = 0.6294 j^^. We see that the scaled data collapse 
into a scaling-function curve. We again confirm that the phase transition belongs to the 
three-dimensional-Ising universality class. In addition to this, from the limiting value of 
the high-temperature side of the scaling function, we estimate the critical amplitude as 
A^"*" = 2.09(13) for k = 2 and j = 0.3; more specifically, we read off the value of = 13 
around the regime {T-T^)L^''' ^ 30, and as for an error indicator, we accepted the amount 
of the data scatter among = 5, . . . , 13. 

Similarly, we determined A^"^ for various parameter ranges of both j and /t. In Fig. 
we plotted the amplitude A^"*" for k = 1, 2 and 4 with A(= j + 0.25) varied. In the plot, we 
observe a clear signature of the power- law singularity as described by Eq. (fT^ . Hence, we 
confirm that the cross-over behavior (fT^ is realized actually around the multi-critical point 
j = —0.25. Moreover, in the figure, we notice that the data for k, = 1,2 and 4 almost overlap 
each other. It would be rather remarkable that the amplitude A^"*" itself hardly depend on the 
parameter k. This fact indicates that the multi-criticality, namely, the singularity exponent 
{—i' + u)/(f), stays universal with respect to the self-avoidance parameter k. Such univers ality 
was first reported by the series-expansion analyses surveying the range ofK = 0.5,...,3 |l7|. 

From the slopes in Fig. we obtained the singularity exponent as (— z> + z/)/0 = 0.422(6), 
0.405(5). and 0.415(7) for k = 1, 2 and 4, respectively. We estimate the singularity exponent 
as (— z> + z/)/0 = 0.415(20) consequently. 

Let us mention some remarks on this estima te ( — z> + ]/)/(/) = 0.415(20). First, this result 
excludes such a possibility z> > z/ as z> = 1.2(1) 15|. Rather, our result supports the results 



of z/ = 0.44(2) {k = 1) with the Monte Car 
the low-temperature-series-expansion result 



o method Q] and u = 0.46(1) (k = 1) with 



17| . (The latter is obtained from a = 0.62(3) 



(k = 1) together with the hyperscaling relation a = 2 — dp.) Note that our preliminary 
survey in the preceding subsection also indicates a signature oi v < v. 

Second, postulating the value z> 0.45 close to the aforementioned existing values, we 
obtain an estimate for the crossover exponent cf) ^ 0.43. The present result contradicts 
the result = 1.1(1) Q| determined with the cluster- variation method. In the succeeding 
section, we will provide further support to ~ 0.43, performing the multi-critical scaling 
analysis based on the relation (fT!?|) . 
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C. Multi-critical scaling analysis 



In the above, we obtained an estimate for the crossover exponent ^ 0.43 from the 
power-law singularity of the amplitude A^"*", accepting the value z> ^ 0.45 advocated by 



Refs. 



17[. In this subsection, we provide further support to these exponents. We carry 



out an multi-critical (crossover) scaling analysis based on Eq. ()13|). For finite size L, the 
scaling-hypothesis formula should be extended to, 

i = LX{{T -T^)L^/\^L'^/^). (14) 

Based on this formula, in Fig. IHl we present the scaled data, (T — Tc)L^^''-^/ L, with fixed 
AL'^/'^ = 2 and k = 2. Here, we set the exponents z> = 0.4 and (p = 0.6 for which we found 
the best data collapse. Surveying the parameter space beside this condition, we obtained the 
critical exponents as z> = 0.45(15) and = 0.6(2). These estimates agree with the analysis 
in the preceding subsection. 

We stress that the use of ^ greatly simplifies the scaling analyses, because ^ has a fixed 
scaling dimension, namely, [length]^. For instance, as for other quantities such as the sus- 
ceptibility, we need to determine the exponent 7 in addition to i>. In this sense, the present 
approach via the transfer matrix is advantageous over other approaches. 



IV. SUMMARY AND DISCUSSIONS 

We investigated the (multi) criticality of the gonihedric model with the extended 
parameter space (j21). The model is notorious for its slow relaxation to the thermal equi- 
ibrium (glassy behavior), which deteriorates the efficiency of the Monte Carlo sampling 
22j |. Aiming to surmount the difficulty, we employed the transfer-matrix method. We im- 
plemented Novotny's idea extending it so as to incorporate the plaquette-type 
interactions (Sec. |n}. The present approach enables us to treat arbitrary number of spins 
per one transfer- matrix slice, admitting systematic finite-size-scaling analyses; see Fig. |3]for 
instance. 

The transfer-matrix calculation has an advantage in that it yields the correlation length 
immediately. Because the correlation length has a known (fixed) scaling dimension, the sub- 
sequent scaling analyses are simplified significantly. Moreover, with the correlation length, 
we are able to calculate the Roomany-Wyld approximate beta function (T) ((TUJ. With 
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use of P§^(T), we surveyed the critical branch of the phase diagram (Fig. Thereby, 
we observed that the criticahty is maintained to be the three-dimensional-Ising universal- 
ity class all along the phase boundary. On closer inspection, we found an indication of 
a crossover critical phenomenon such that the slope of /5^^(T) in the off-critical regime, 
typically, T — > 3 (j = —0.05), acquires a notable enhancement. This fact indicates 
that a multi-criticality with smaller z> emerges as j —0.25. This observation supports the 
claim [l^ that the nonstandard criticality reported so far [l5, 16, 17, 18, 19, 20 1 could be 
attributed to the end-point criticality specific to j = —0.25. 

Aiming to clarify the nature of this multi-criticality, we analyzed the end-point singularity 
of the amplitude of the correlation length (jl2|) . As shown in Fig. El the amplitude 
exhibits a clear power-law singularity, from which we obtained an estimate for the singularity 
exponent (— z> + = 0.415(20). This result supports the above-mentioned observation 



that an inequality v < v shou 
oi v > V advocated in Ref. [15 



d hold, and in other words, it excludes such a possibility 
]. Rather, our result supports the Monte Carlo simulation 
result z> = 0^.44(2) (k = 1) [iq and the low-temperature-series-expansion result z> = 0.46(1) 
(k = 1) |17]. Postulating z> ^ 0.45, we arrive at an estimate for the crossover exponent 
(j) 0.43. This exponent is to be compared with the result (f) = 1. 1(1) determined with the 
cluster- variation method The discrepancy between the result ^ and ours seems to be 
rather conspicuous. 

We then carried out the multi-critical scaling analysis (|Tl|l in order to provide further 
support to our estimate ~ 0.43 based on z> ^ 0.45. We found that a good data collapse 
is attained for z> = 0.4 and (f) = 0.6 under k = 2 and AL'^^'^ = 2 (Fig. IH}. Surveying the 
parameter space, we obtained the estimates = 0.6(2) and z> = 0.45(15). These exponents 
agree with the above-mentioned analysis via the critical amplitude A^"*". 

As a consequence, we confirm that the whole analyses managed in this paper lead a self- 
consistent conclusion. Regarding the discrepancy on 0, we suspect that the value = 1.1(1) 
[ifll might be rather inconceivable. Nevertheless, in order to fix the multi-criticality more 
definitely, further elaborate investigations would be required. As a matter of fact, a possible 
slight deviation of the multi-critical point from j = —0.25 was ignored throughout the 
present work as in Ref. Justification of such a treatment might be desirable. In any 

case, the present approach, which is completely free from the slow-relaxation problem, would 
provide a promising candidate for a first-principles-simulation scheme in future research. 
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FIG. 1: (a) A schematic phase diagram for the gonihedric model is shown For large k, A 
second-order phase transition occurs. The criticality has been arousing much attention, (b) For an 
extended parameter space, Eq. 0, there emerge rich phases accompanying a multi-critical point 
[2^ : here, the self-avoidance parameter k is fixed {k = 1). In terms of this extended parameter 
space, the transition point in Fig. ^ (a) is identified with the multi-critical point at j = —0.25. 
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FIG. 2: Novotny invented a new scheme to construct the transfer matrix (TM) 
which allows us to treat arbitrary number of spins per one TM slice. We extend his scheme to 
incorporate the plaquette-type (next-nearest-neighbor and four-spin) interactions, aiming to treat 
the gonihedric model The contributions from the "leg," "planar," and "rung" interactions are 
considered separately; see Eq. (jSJ. With use of the translation operator , we build a bridge 
between the \/iV-th neighbor spins along the leg (inter-leg interaction). 
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FIG. 3: The beta function (3^^ (T) ^ is plotted for k = 2 and various j. For a comparison, we 
presented a slope (dashed line) corresponding to the three-dimensional-Ising universality class (z^ = 
0.6294 [3^); we see that the criticality is maintained to be the three-dimensional-Ising universality 
class for a wide range of j. In fact, from the slopes at the fixed points of /3^^(T), we obtain 
an estimate for the correlation-length critical exponent u = 0.638(5); see text for details. Inset: 
Plotting the zero-points of (3i^ (T), we determine a phase boundary Tc{j), which corresponds to 
the critical branch in Fig. ^(b). 
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FIG. 4: Scaling plot for the correlation length, namely, {T — Tc)L^^'^-(,\T — T^'^ , is shown for k = 2 



36|. 



and j = 0.3. Here, we postulated the three-dimensional-Ising universality class u = 0.6294 
The symbols, +, x, *, □, and ■, denote the system sizes of = 5, 7, 9, 11, and 13, respectively. We 
confirm that the transition belongs to the three-dimensional-Ising universality class. Furthermore, 
from the plateau in the high-temperature side, we obtain an estimate for the critical amplitude 
N+ = 2.09(13); see text for details. 
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FIG. 5: Correlation-length critical amplitude A'^"'" ()12|) is plotted for various A(= j + 0.25) and 
K = 0.5, 1 and 4. The symbols, +, x, and *, stand for the self-avoidance parameter k = 0.5, 1, 
and 4, respectively. The data indicate a clear power-law singularity, Eq. ()12[). From the slopes, we 
estimate the singularity exponent as (— z/ + i')/(p = 0.415(20). 
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FIG. 6: Multi-critical (crossover) scaling plot (HH), (T - Tc)L^/''-^/ L, for k = 2 and AL'^/'^ = 2 is 
shown. Here, we set u = 0.4 and (j) = 0.6, for which we found the best data collapse. The symbols, 
+, X, *, □, and ■, denote the system sizes of = 5, 7, 9, 11, and 13, respectively. 
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